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Assessing the correspondence between spontaneous and stimulus-driven neural activity can reveal intrinsic 
properties of the brain. Recent studies have demonstrated that many large-scale functional networks have a 
similar spatial structure during spontaneous and stimulus-driven states. However, it is unknown whether 
the temporal dynamics of network activity are also similar across these states. Here we demonstrate that, in 
the human brain, interhemispheric coupling of somatosensory regions is preferentially synchronized in the 
high beta frequency band (—20-30 Hz) in response to somatosensory stimulation and interhemispheric 
coupling of auditory cortices is preferentially synchronized in the alpha frequency band (—7-12 Hz) in 
response to auditory stimulation. Critically, these stimulus-driven synchronization frequencies were also 
selective to these interregional interactions during spontaneous activity. This similarity between 
stimulus-driven and spontaneous states suggests that frequency- specific oscillatory dynamics are intrinsic 
to the interactions between the nodes of these brain networks. 

Neurons, in vivo, are active in the absence of external input. Increasing evidence suggests that this 
spontaneous activity has a coherent structure and may play a role in a host of cognitive and neural 
processes, such as motor performance 1 , learning and memory 2 , and neural development 3 . In support of 
this possibility, recent neurophysiological studies have shown that the activity of individual neurons is similar 
across spontaneous and stimulus-driven brain states 4 " and a recent transcranial magnetic stimulation study 
suggested the same 9 . This similar activity has motivated a number of hypotheses about the mechanism by which 
spontaneous neural activity contributes to cognition 4-6,9-11 . 

The neurophysiological evidence from single neurons has been complemented by findings from neuroimaging 
studies that have revealed a correspondence between spontaneous and stimulus-driven states regarding the 
spatial structure of large-scale functional networks, consistent with the idea that the brain has an intrinsic 
functional architecture 12-15 . This spatial correspondence has been based on interregional correlations of spon- 
taneous, slow fluctuations in the brain's hemodynamic activity (< .1 Hz) 12-15 with some evidence that these 
fluctuations may be associated with gamma band electrophysiological activity (30+ Hz) 16-18 or that the precise 
frequency associated with these slow fluctuations differ across the brain 19-20 . However, it is unknown whether the 
frequency bands associated with interregional synchrony during stimulus-driven and spontaneous states of these 
large-scale interactions are related 9 . 

Current evidence suggests that stimulus -driven interregional synchrony is associated with temporal dynamics 
in a broad range of frequency bands (here termed "oscillatory synchronization") that may vary as a function of 
task and brain network. For example, spatial navigation tasks modulate hippocampal neuronal synchrony in the 
theta frequency band (4-8 Hz) 21 , and median nerve stimulation modulates neural synchrony in the somatosen- 
sory cortex in the beta frequency band (13-30 Hz) 22,23 . Here we test the hypothesis that the frequency bands 
associated with stimulus-driven oscillatory synchronization between specific neural regions are also characteristic 
of their spontaneous synchrony. 

Results 

We used magnetoencephalography (MEG) to compare stimulus-driven and spontaneous oscillatory phase lock- 
ing 22,24,25 - a measure of the precise temporal synchronization of neural signals - in somatosensory and auditory 
regions. The regions chosen for phase locking analysis were in different hemispheres of the brain and further than 
6 cm apart to avoid potential artifactual crosstalk that occurs between relatively proximal locations in MEG 24 . 
Data for evaluating stimulus-driven oscillatory synchrony in the somatosensory regions were obtained by having 
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Somatosensory and motor event related field 




B Auditory event related field 




Figure 1 | (A) Single subject example of the event related fields following a 
button press with the left index finger (right hemisphere image) or the 
right index finger (left hemisphere image). Noise normalized, dynamical 
statistical parametric mapping images are plotted 55 . The regions chosen 
for LM1, SI, and S2 are labeled. (B) Single subject example of the event 
related fields during passive listening to auditory clicks. The regions chosen 
for LAud and RAud are labeled. 



subjects (N = 10) perform a cued bimanual button press task (see 
figure 1A for an example event related response to button press). 
Phase locking was calculated between a region in the right hemi- 
sphere primary somatosensory cortex (RSI) and the left hemisphere 
secondary somatosensory cortex (LS2) identified via analysis of the 
data from 0-150 ms after the onset of the button press (figure 1A; 
synchrony between RSI and LSI could not be reliably evaluated 
because they are less than 6 cm apart). Data for evaluating stimu- 
lus-driven oscillatory synchrony in the auditory cortices were 
obtained while subjects passively listened to clicks (15, 20, 30, or 
40 Hz, 500 ms duration, interstimulus interval =1.5 sec) delivered 
binaurally (see figure IB for an example event related response to 
auditory clicks). Phase locking was calculated between left and right 
hemisphere primary auditory cortices (LAud and RAud), identified 
via analysis of the passive listening data averaged across all frequency 
conditions (figure IB). Data for evaluating spontaneous oscillatory 
synchrony were obtained by having subjects fixate a centrally located 
cross for 12 minutes prior to their participation in the button press or 
passive listening task. Spontaneous phase locking was evaluated 
between the aforementioned interhemispheric somatosensory and 
auditory sites during the passive fixation period. All statistical values 
were based on a permutation test and corrected for multiple fre- 
quency comparisons 26 . 



Stimulus-driven and spontaneous synchrony between RSI and 
LS2. In somatosensory cortices, we found increased interhemis- 
pheric phase locking between RSI and LS2 that was restricted to 
the high beta frequency range during the button press relative to a 
prestimulus baseline (figure 2A; p = .008 from 20-28 Hz). This 
finding is consistent with previous studies that demonstrated that 
phase locking increases in the beta frequency band between RSI and 
LS2 following median nerve stimulation 22,23 and reports indicating 
the importance of beta band synchrony in the sensory- motor sys- 
tem 27 ~ 2S . The similarity between our phase locking results to previous 
reports of synchrony during somatosensory stimulation 22,23 and the 
time window of analysis (0-150 ms post button press onset) suggest 
that this RS 1 -LS2 phase locking reported is a result of proprioception 
during the button press. Importantly, significantly greater phase 
locking in the high beta frequency range was also seen for spon- 
taneous activity between these same sites (RSI and LS2), relative to 
both noise and the RH motor cortex (RM1)-LS2 control site pair 
(figure 2B;p< .002 from 22-30 Hz). The largest peak in the RSI -LS2 
phase locking occurred in the alpha frequency band, though the 
synchrony at these frequencies was not significantly different be- 
tween RS1-LS2 and RM1-LS2. Indeed, the somatomotor rhythms 
have both alpha and beta frequency components 30 and this alpha 
synchrony may be involved in somatosensory-motor interactions. 
It is noteworthy that a peak of phase locking in the high beta fre- 
quency band is also visible during the prestimulus period (figure 2A), 
consistent with the idea that this frequency is important to spon- 
taneous synchrony between RSI and LS2. 

Stimulus-driven and spontaneous synchrony between RAud and 
LAud. A similar pattern of correspondence between the stimulus- 
driven and spontaneous interhemispheric oscillatory synchrony was 
found between the auditory sites. Specifically, we found increased 
phase locking between RAud and LAud that was restricted to the 
alpha frequency range during the presentation of the clicks relative to 
a prestimulus baseline (figure 3 A; p < .05 from 7-15 Hz). While we 
are unaware of any previous studies that examined phase locking 
between RAud and LAud in response to auditory stimulation, 
previous work has shown that the local neural response in the 
alpha frequency band in both RAud and LAud is modulated by 
auditory stimuli 3133 . These same sites (RAud-LAud) also showed 
significantly greater phase locking in the alpha frequency range for 
spontaneous activity, relative to both noise and the RAud-LMl 
control site pair (figure 3B; p = .006 from 10-14 Hz). The fact 
that a peak of phase locking in the alpha frequency band is also 
visible during the prestimulus period (figure 3A), further supports 
the idea that this frequency is important to spontaneous synchrony 
between RAud and LAud. 




Frequency (Hz) Frequency (Hz) 

Figure 2 | Stimulus-driven and spontaneous synchronization between left and right hemisphere somatosensory cortices. (A) Phase locking between 
RS 1 and LS2 with respect to frequency during a cued button-press task. The phase locking during the button-press and during the period prior to 
presentation of the button-press cue (pre-stimulus baseline) is shown. Significantly greater phase locking was seen during the button-press compared to 
the pre-stimulus baseline period only in the high beta frequency range (20-28 Hz indicated by the grey bar above the data). (B) Spontaneous phase 
locking between RSI and LS2 and between RM1 and LS2 (control pair) with respect to frequency. RM1-LS2 were chosen as the control pair because their 
Euclidean distance is the same as between the RS 1 -LS2 pair of interest. Significantly greater phase locking was seen for RS 1 -LS2 versus noise and RM 1 -LS2 
control pair regions only in the high beta frequency range (22-30 Hz indicated by the grey bar above the data). 
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Frequency (Hz) Frequency (Hz) 

Figure 3 | Stimulus-driven and spontaneous synchronization between left and right hemisphere auditory cortices. (A) Phase locking between RAud 
and LAud with respect to frequency during presentation of auditory clicks and during the period prior to presentation of the clicks (pre-stimulus 
baseline). Significantly greater phase locking was seen during auditory stimulation compared to the pre-stimulus baseline period only in the alpha 
frequency range (7-15 Hz indicated by the grey bar above the data). (B) Spontaneous phase locking between RAud and LAud versus between RAud and 
LM1 (control pair regions) with respect to frequency. RAud-LMl were chosen as the control pair because their Euclidean distance is the less than the 
RS1-LS2 pair of interest. Significantly greater phase locking was seen for RAud-LAud than for RAud-LMl only in the alpha frequency range 
(10-14 Hz indicated by the grey bar above the data). 



Slow fluctuations of spontaneous synchrony. These findings estab- 
lish that the frequencies associated with enhanced, interhemispheric 
neural synchrony between somatosensory regions due to proprio- 
ception during a motor task, and between auditory cortices when 
listening to clicks, are also critical to spontaneous synchrony within 
these systems. To address whether our findings might be related to 
the slow fluctuations (< .1 Hz) of correlated activity seen in neuro- 
imaging studies of spontaneous functional connectivity 12,13,15 , we 
examined the rate at which this relatively high frequency inter- 
regional synchronization varied over time. Specifically, we used a 
Fourier analysis to examine the rate of change over time of the 
phase locking at the target frequency of both interhemispheric 
interactions (figures 4A and 4C). For the two pairs of regions, RS1- 
LS2 and RAud-LAud, a greater amount of the change over time of the 



high beta and alpha phase locking, respectively, was accounted for by 
lower frequencies than higher frequencies, with the lowest 
frequencies (e.g. < .1 Hz) accounting for the most variability 
(RS1-LS2: mean R = -.50, t = -5.62, p < .001; figure 4B; RAud- 
LAud: mean R = -.60, t = -13.44, p < .001; figure 4D). Thus, 
spontaneous interregional synchronization between these brain 
areas is characterized by very slow fluctuations of higher frequency 
oscillations 19,20,34 . Critically, these higher frequency oscillations are 
the same for stimulus-driven and spontaneous synchrony. 

One question is whether these slow fluctuations of higher fre- 
quency phase locking are primarily driven by the regions slowly 
coupling and decoupling or whether they are governed by increases 
and decreases in local activity. While increases and decreases in local 
activity do not directly affect phase locking, they modulate the 
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Figure 4 | (A) Spontaneous time courses of the 24 Hz phase locking between RSI and LS2 for a representative subject. As illustrated, synchrony 
(phase-locking values) fluctuated very slowly over time (on the order of sees). Phase locking is plotted relative to the average of 1000 noise simulations 
using random, Gaussian-distributed noise. (B) The Fourier spectrum of the 24 Hz synchrony between RSI and LS2. Power increased linearly from the 
highest to the lowest frequencies. A significant (p < .01) negative linear relationship between frequency and power was seen in 8/10 subjects with an 
overall highly-significant relationship across subjects (mean R = —.50, t = —5.62, p < .001). Furthermore, the slope of the 24 Hz spectrum was 
significantly steeper (more negative) than random noise (t = —2.14, p =.03) while the 18 Hz spectrum (not pictured) did not differ from noise 
(t = — 1.31, p > .1) indicating that the high beta frequency synchrony between RSI and LS2 fluctuates very slowly. (C) Spontaneous time courses of the 
10 Hz phase locking between RAud and LAud for a representative subject illustrating that the synchrony varied slowly over time. (D) The Fourier 
spectrum of the 10 Hz synchrony between RAud and LAud. As in panel B, power increased linearly from the highest to the lowest frequencies, and there 
was a significant negative linear relationship between frequency and power in 9/10 subjects with an overall highly-significant relationship across subjects 
(meanR= — .60,t= —13.44, p < .001). Also, the slope ofthe 10 Hz spectrum was significantly steeper than random noise (t = — 3.03, p =.007) while the 
20 Hz spectrum (not pictured) was not (t = — 0.99, p > . 1 ) indicating that the alpha frequency synchrony between RAud and LAud changes very slowly. 
Note that while the data B and D were well fit by a negative linear relationship, they were also well fit by a 1/frequency relationship. 
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signal-to-noise ratio, which does affect measures of phase locking 35,36 . 
To address this question we examined the correlation over time 
between local power and phase locking. A significant correlation 
was seen between RSI power and RS1-LS2 phase locking at 24 Hz 
(mean R = .17, t = 8.32, p < .001), LS2 power and RS1-LS2 phase 
locking at 24 Hz (mean R = .13, t = 3.23, p < .05), RAud power and 
RAud-LAud phase locking at 10 Hz (mean R = .12, t = 3.14, p < 
.05), and LAud power and RAud-LAud phase locking at 10 Hz 
(mean R = .20, t = 16.4, p < .001). Additionally, the RSI and LS2 
power-phase locking correlations were significantly greater at 24 Hz 
thanat 18 Hz(p < .01 in both cases) and the RAud and LAud power- 
phase locking correlations were significantly greater at 10 Hz than at 
20 Hz (p < .005 in both cases). Finally, it should be noted that the 
power correlation between RSI and LS2 at 24 Hz was not signifi- 
cantly larger than at 18 Hz (t = 1.91, p =.09) nor was the power 
correlation between RAud and LAud at 10 Hz significantly larger 
than at 20 Hz (t = 2.03, p =.07), though both showed a trend 
towards significance. These results demonstrate that some of the 
fluctuations over time of the interregional phase locking can be 
explained by increases and decreases in the local power in the 
respective areas. However, the correlation coefficients do not exceed 
.20, indicating that these local changes do not explain a substantial 
proportion of the variance (< 5%) of the phase locking over time. 
Thus, our results suggest that RS1-LS2 and RAud-LAud slowly go in 
and out of synchrony during the spontaneous state at 24 Hz and 
10 Hz respectively, at least partially independent of local activity 
changes. 

Discussion 

Our results suggest that activity in the same frequency bands under- 
lies interregional synchrony for both stimulus -driven and spontan- 
eous brain states. Specifically, increased interhemispheric phase 
locking in the beta frequency band between RSI and LS2 and in 
the alpha frequency band between RAud and LAud was associated 
with somatosensory and auditory response respectively. Further- 
more, interhemispheric phase locking in these same frequency bands 
was characteristic of spontaneous communication between these 
regions. This correspondence between stimulus-driven and spontan- 
eous oscillatory synchronization suggests that the dynamics are not 
purely state dependent, appearing and disappearing as needed. 
Rather, our results suggest that the dynamics at these frequencies 
are intrinsic to these interregional interactions. 

In support of the hypothesis that frequency specific synchrony is 
intrinsic to interregional interactions, physiology and neural model- 
ing suggests that relatively stable properties of neural regions, such as 
different conduction delays between brain regions and other network 
properties can guide oscillatory dynamics 37 40 . One putative role for 
these intrinsic dynamics is to allow for dynamic control over how 
and where information is gated 38,41 . Further studies will be required 
to determine whether frequency specific oscillatory dynamics extend 
beyond interhemispheric coupling of early sensory cortices. 

The present results suggest that interacting brain regions sponta- 
neously synchronize in particular frequency bands and that stimu- 
lus-driven activity increases the interregional synchrony at that 
frequency. These intrinsic oscillatory dynamics are generally consist- 
ent with the hypothesis that stimulus-driven activity can be described 
as a perturbation of the ongoing, spontaneous electrophysiological 
brain activity 11,42 . Most previous studies have examined these pertur- 
bations relative to the prestimulus period (i.e. in the context of the 
task or stimulus). The similarity between the prestimulus baseline 
and the spontaneous activity seen in the present results extend this 
framework beyond a task or stimulus context to the resting-state. 
Furthermore, most previous studies have examined these perturba- 
tions in local neural activity 11,43 . Our results suggest that a similar 
mechanism may underlie interregional synchrony at the intrinsic 
oscillatory frequencies described here. 



This is not meant to imply that all types of oscillatory synchrony 
are intrinsic to particular brain interactions. For example, gamma 
frequency dynamics are seen ubiquitously throughout the brain. 
Gamma synchrony is associated with distinct cognitive states such 
as attention 44 , learning 45 , and working memory 46 , processes that 
inherently involve many brain networks. Thus, networks may be 
partly characterized by intrinsic oscillatory mechanisms while also 
being subject to transient or global oscillatory states. Indeed, 
dynamics in multiple oscillatory frequencies can exist within net- 
works and often interact with one another, leading to frequency 
nesting and cross-frequency interactions 47,48 . 

It is unclear how the present results relate to previous reports 
emphasizing the role of gamma band oscillations in slow spontan- 
eous activity 1618 because these previous studies did not directly com- 
pare stimulus-driven and spontaneous synchrony as was done here. 
The present results demonstrate that the spontaneous state is char- 
acterized by slow changes in higher frequency synchrony, and that 
these higher frequency dynamics are similar in stimulus-driven and 
spontaneous states. At long time-scales, this result is consistent with 
neuroimaging studies showing that spontaneous activity is charac- 
terized by slow fluctuations of correlated activity 12,13,15 . At relatively 
brief timescales (on the order of hundreds of milliseconds), this result 
is also consistent with neurophysiological evidence in animals that 
show a marked similarity between the temporal structure of spon- 
taneous and stimulus-driven activity at the spatial scale of individual 
neurons and at the temporal scale of milliseconds 4 8 . Accordingly, 
the present work bridges the gap between these previously disparate 
human neuroimaging and animal neurophysiology studies. 

Although our data do not address the function of ongoing neural 
synchrony, one intriguing possibility is that it provides stability to the 
interregional interactions that underlie brain networks 13 . Consistent 
with this idea, recent modeling work demonstrates that interregional 
synchronization can improve a network's resilience against the 
potentially destabilizing effects of noise 49 . Therefore, sustained oscil- 
latory synchrony may be a way for nodes of a neural network to resist 
noise in order to maintain coherence and stabilize the information 
represented in that network. 

The present data demonstrate that spontaneous interregional syn- 
chrony slowly varies over time (i.e. non-stationarity). These results 
are generally consistent with previous fMRI 50,51 , MEG 19,20 , and inter- 
cranial EEG 16 " 18 studies showing similar effects. These previous 
results have looked primarily at whether interregional correlations 
between the overall level of activity fluctuated over time (though 51 
examined slow changes in phase coupling using fMRI as well). Our 
results demonstrate that the relative phase of the regions also fluc- 
tuates slowly over time, with the phase of RS1-LS2 and RAud-LAud 
slowly coupling and decoupling over the spontaneous period. 
Furthermore, these changes in phase locking were at least partially 
independent of local changes in the power of the activity within these 
regions. 

One potential limitation of these data is that the RS 1 -LS2 stimulus- 
related synchrony relies on somatosensory proprioception during a 
motor task. This task may confound motor and somatosensory activ- 
ity and is not well matched to the passive auditory task used to 
examine RAud-LAud synchrony. However, previous studies using 
passive median nerve stimulation have reported increased synchrony 
specific to the beta-band in these tasks 22,23 , which is consistent with 
the present results despite the very different tasks. Additionally, it is 
not clear if the stimulus-related increases in synchrony between RS 1 - 
LS2 and RAud-LAud were due primarily to increased coupling 
between these regions or due to local changes in activity that result 
from the stimuli 35,36 . Neither of these two possibilities would obviate 
the importance of the relevant frequencies to interregional commun- 
ication, but they would affect the interpretation of how changes in the 
phase locking at these frequencies result from stimulus processing. In 
order to definitively resolve this ambiguity between interregional 
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coupling and local power, new methods will need to be developed 
and tested to assess phase locking independent of local power 
changes. Finally, we only compared stimulus-driven and spontan- 
eous synchrony in two, relatively low-level sensory systems. Further 
studies will be required to determine whether the frequency corres- 
pondence between spontaneous and stimulus driven activity is a 
general property of the brain, particularly for higher-level networks. 

Our findings provide evidence that, at the level of large-scale 
neural interactions, spontaneous neural synchrony has a coherent 
temporal structure, with the nodes of different networks showing 
preferential biases for interacting within limited frequency bands. 
They further suggest that frequency-specific oscillatory dynamics 
are intrinsic to the interactions between the nodes of these brain 
networks. These dynamics can potentially be tied to specific bio- 
logical and functional properties of neurons and networks 13,37 " 40,49 . 
Therefore, a better understanding of intrinsic dynamics and how 
they are altered in pathological conditions may yield critical insights 
into the mechanisms behind spontaneous oscillatory synchrony and 
abnormalities of functional connectivity. 

Methods 

Subjects. A total of 15 individuals participated in the experiments. One subject was 
excluded due to unusually large cardiac and respiratory artifacts and another subject 
was excluded due to head movement in excess of .5 cm, thus 13 subjects were 
included in the analysis (8 males, mean age = 27.3, SD = 8.1). Of these 13 subjects, 
seven individuals participated in both of the stimulus -driven experiments and six 
individuals participated in just one of the two stimulus-driven experiments (a total of 
10 subjects in each experiment). Spontaneous activity was recorded from all subjects. 
All subjects were naive to the goals of the study. The Institutional Review Board of the 
National Institutes of Health approved all procedures and written informed consent 
was obtained for all subjects. 

Recording. Neuromagnetic responses were recorded at 600 Hz using a 275 channel 
whole head MEG system in a shielded room (VSM MedTech, Ltd., Canada). The 
magnetometer is equipped with 275 radial gradiometers (273 were functional) and 
synthetic 3 rd order gradient noise cancellation was used. Head position coils were 
placed at the nasion and left and right preauricular points to coregister the anatomical 
MRI and the MEG sensors. Head position was determined at the beginning and end of 
each run to ensure that head movements did not exceed .5 cm for any subject. 
Eyeblinks and eye movements were recorded using a bi-polar EOG electrode placed 
about each subject's left eye. 

Structural MRI. Structural MRI images were obtained separately using a 3-Tesla 
whole-body scanner {GE Signa, USA). A high -resolution Tl-weighted 3D volume 
was obtained for each subject. Freesurfer™ was used to create a cortical surface model 
for each subject using an automatic reconstruction algorithm. The cortical white 
matter was segmented providing a topologically correct representation of the surface 
with approximately 150,000 vertices per hemisphere. The cortical surface was then 
decimated to approximately 4,000 source dipoles per hemisphere, approximately 1 
dipole every 10 mm along the cortical surface. The MEG data were coregistered with 
anatomical MR images using fiduciary headpoints. 

Behavioral tasks. Each subject first participated in a 12 minute scan of spontaneous 
activity where their task was to fixate on a centrally presented cross. Subjects were not 
given any details about later tasks during this scan. Following the recording of 
spontaneous activity, the subjects were presented with a series of button-press trials. 
Every 2 seconds, a white fixation cross in the middle of the screen would change to 
green or purple for 200 ms. The subjects were instructed to press a button with their 
left index finger when the cross turned green and a different button with their right 
index finger when the cross turned purple. The subjects were presented with 100 trials 
in each condition. Following this task, the subjects were presented with a series of 
auditory clicks trials binaurally using MEG compatible ear buds. Each trial consisted 
of 500 ms of a 1910 Hz pure tone modulated by square waves at 15, 20, 30, or 40 Hz 
followed by 1500 ms of silence. The subjects were presented with 100 trials of each 
frequency and asked only to fixate on a centrally presented cross throughout the 
experiment. 

MEG analysis. Artifact removal. For artifact removal, we first visually inspected the 
data for any respiratory artifacts. These artifacts are generally uncommon in MEG 
data, particularly with the use of third-order gradiometer compensation. However, if 
subjects accidentally have metal on or in their persons these artifacts can be 
substantial (one subject in our study was excluded for this reason). We then used a 
short-time Fast Fourier filter to band-pass the data between 1 and 50 Hz (the filter 
uses the overlap-add technique with a time window of 4096, the highpass transition 
width is 0.0 Hz and the lowpass transition width is 5.0 Hz). This removed low 
frequency drift, any residual respiratory artifact, 60 Hz line noise, and any DC offset. 



Cardiac artifacts were removed using an independent component analysis -based 
procedure following Liu et al. 52 . Briefly, the MEG sensor data were decomposed into 
273 independent components (ICs) using EEGLAB 53 . ICs were identified as being a 
cardiac artifact if the IC had a peak in its autocorrelation spectrum between .6-1.5 Hz 
(by observing whether there were regular peaks in the autocorrelation every 667- 
1667 ms), the IC had a timecourse that contained periodic features that were similar 
to those seen on an electrocardiogram, and had power in the MEG sensors that 
typically contain cardiac features. Across our subjects, between 1 and 3 cardiac ICs 
were rejected and the remaining ICs were reassembled for further processing. 

To monitor ocular muscle activity we measured EOG along with the MEG mea- 
surements. The EOG timecourses were pseudo-Z transformed into standard devi- 
ation (SD) units (by subtracting the mean across the epoch and dividing by the SD). 
These data were visually inspected and the minimum size of each subject's eye blinks 
in these units was determined (mean — 1.67 [in SD units], min 1.5, max 2.0). 
Centered about each point that exceeded these thresholds, the data 300 ms before and 
500 ms after these points were excluded from further analysis (on average this 
removed 1200 ms per eye blink). These additional windows of data were excluded so 
that sufficient pre- and post-artifact data were removed to ensure no artifacts 
remained in the analyzed timecourses. These conservative thresholds and windows 
removed on average 38% of the time points across subjects, leaving 7.4 minutes of 
artifact-free data on average. Note that simulations demonstrate that the discon- 
tinuities in the data created by this artifact removal procedure do not substantially 
affect the calculated phase locking above 3 Hz 24 . 

Inverse solution. The exact location of the cortical current sources cannot be precisely 
determined using the measured magnetic fields from outside the head, and therefore 
it is estimated using cortically constrained minimum norm estimate (MNE), 
described extensively elsewhere 54 ' 55 . Briefly, a linear inverse operator W is applied to 
the measured signal to calculate the MNE 

y(t)=Wx(t) 

where x(t) represents the MEG channel data at time t and y(t) is the corresponding 
current projected onto the cortical surface. The expression of W is defined as 

W — RA T (ARA T + A Z C) ~~ 1 . 

where C and R are the noise and source covariance matrices respectively. A defines the 
transformation from orthogonal unit current dipoles to measured magnetic fields and 
uses the boundary element method 56 57 . A, is a weighting factor that is used to avoid the 
magnification of errors in the data and X 2 ~ (1/SNR). We used a value of 3 for this, as 
is typically done in MEG analysis (Hamalainen, MNE software user's guide version 
2.7, 2009). Furthermore, because cortical neurons are known to be preferentially 
oriented perpendicular to the cortical surface, we used a loose orientation constraint. 
Specifically, the component of R normal to the surface was multiplied by 1 and the 
components transverse to the surface was multiplied by .4. To compensate for the bias 
towards superficial currents of the MNE a scaling factor (i.e. depth weighting) of .8 is 
applied to R 22 . The noise covariance matrix was calculated from 12 minutes of con- 
tinuous empty room MEG measurements collected immediately prior to putting the 
subject in the scanner 24 . 

To estimate the time course and statistical significance of the cortical MEG activity, 
noise normalized values were calculated at each time point and each dipole location 55 . 
This transforms power values into dynamic statistical parametric values (dSPM) and 
makes the point spread function of the estimated signal relatively uniform across 
cortical dipoles 58 . For ROI localization these dSPM values are used to describe the 
neural activity. 

ROI localization. To localize LH primary motor cortex (LM1) and LH secondary 
somatosensory cortex (LS2), the locus of peak activity about the central sulcus 
occurring for activity averaged from 0 to 150 ms after subjects pressed the button with 
their left hand was calculated. This time period was chosen because the neural 
response at this time includes somatosensory activity in response to proprioception of 
the button press. To localize LM1 and LS2, the locus of peak activity about the central 
sulcus occurring for activity averaged from 0 to 150 ms after subjects pressed the 
button with their right hand was calculated. The portion of this activity anterior to the 
central sulcus, but posterior to the superior frontal gyrus and middle frontal gyrus, 
was used as the location of LM1 and the portion of the activity just ventral to the 
central sulcus extending into the dorsal insula, was used as the location of LS2 for each 
subject. This time period was chosen to correspond approximately to the 
proprioception in response to the button press. The portion of this activity anterior to 
the central sulcus was used as the location of the RM1 and the portion of this activity 
posterior to the central sulcus was used as the location of RSI for each subject. 
Synchrony between RSI and LSI or between RSI and RS2 could not be reliably 
evaluated due to the potential for artifactual cross-talk between locations less than 
6 cm apart in MEG 24 . 

To localize LH and RH auditory cortex (LAud and RAud), the activity between 50 
and 550 ms after the onset of the auditory clicks was averaged. The peak of the activity 
within a 2 cm radius circle on the surface around Heschl's gyrus (determined based 
on anatomical landmarks provided by Freesurfer™) was used as the location of the 
LAud and RAud ROI for each subject. 

Stimulus-driven phase locking analysis. To determine the trial- by- trial phase 
locking between neural regions, we employed the dynamic spectral statistical 
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parametric mapping method 22 , a method to measure phase synchrony between 
signals projected onto the cortical surface. This method employs the anatomically 
constrained MNE inverse solution 54,55 to determine phase locking values {PLVs) 25 
between regions on the cortical surface. 

To calculate PLVs, the MEG sensor data were filtered using a continuous Morlet 
wavelet transform at each frequency of interest described by the equation: 

G(f/)= ^/- 



where c is the SD of the Gaussian envelope of the wavelet in the time domain. To 

7 

ensure stability of the wavelet transform here we set a to be — - . The wavelet 

2 n/ 

representation of each trial was then mapped from the sensors onto the cortex using 
the MNE inverse solution. The phase was then extracted from the wavelet data 
averaged across all dipoles in each region of interest for each trial about each time 
point and at each frequency of interest. Prior to averaging, the phases were rectified to 
correct for the sign inversion that occurs within large ROIs for due to crossing a sulcus 
or gyrus. Specifically, we iteratively flipped the sign of the dipoles in an ROI until no 
pair of dipoles showed strong negative correlations (r < —.3). The PLV was then 
determined using the following formula: 



PLV(t) - ; 



1 



^expO(<Mt,n)-«> 2 (t,n)) 



where (j>( t,n) is the phase of the signals from the two ROIs at time t and trial n. PLVs 
range from 0 to 1, where 1 indicates perfect synchrony across trials and 0 indicates the 
phase relationship is completely random. The same timeframes used to select the 
ROIs from the evoked signal were used for the stimulus -driven PLVs and equal time 
windows from the prestim were used for prestimulus PLVs. Stimulus -driven PLVs 
were scaled by the PLV expected by chance calculated using the Rayleigh distribution. 

Spontaneous phase locking. Spontaneous phase locking measures the variability 
over time of the phase difference between the seed and the other cortical locations 24 . 
Specifically, the spontaneous PLV is defined as: 



PLV i = 



N 



where N is the number of time points in the sample and 6 see d(n) and G](n) are the 
phase of the wavelet convolved data in the seed ROI and the cortical locations 'T 
respectively. Spontaneous PLVs were scaled by the PLV found in empty room noise as 
described in Ghuman et al. 201 1 24 . 

The Fourier spectrum of spontaneous phase locking. To assess the rate of change of 
the spontaneous phase locking, we first calculated the phase locking in 1 second, half- 
overlapping windows {i.e. the N in the equation above was 1 second for each window) 
to across the artifact free portion of the 12 minute spontaneous activity scan 
(e.g. figures 4a and 4c). This yielded a timecourse of the PLVs over the spontaneous 
activity scan. We then calculated the Fourier spectrum of this PLV timecourse in 
Matlab™. 24 Hz was chosen for the RS1-LS2 interaction {18 Hz was chosen as a 
control frequency) and 1 0 Hz was chosen for the RAud-LAud interaction (20 Hz was 
chosen as a control frequency) because these were the frequencies that showed the 
greatest spontaneous phase locking relative to baseline (see figures 2b and 3b). The 
procedure above was then performed on 1000 simulations of Gaussian distributed 
random data to estimate the time- windowed PLVs and Fourier spectra in noise. The 
fit of the Fourier spectrum of the resting-state and noise PLV timecourses to a line and 
a 1/f curve using a least squares fit was determined. 

In addition to calculating the phase locking in these windows, we calculated the 
local power in RSI andLS2 at 24 and 18 Hz, and RAudand LAudat 10 and 20 Hz by 
taking the absolute value of the wavelet transformed response in each region. We then 
calculated the correlation between these local power time courses and the phase 
locking time courses described in the previous paragraph. This allowed us to examine 
the potential effect of local power changes over time on the interregional phase 
locking 3536 . 

Statistics. Nonparametric statistics were then used to compare conditions and 
control for multiple comparisons 26 . All the frequency points that were p < .05 for the 
prestimulus vs. the poststimulus periods were determined and clustered on the basis 
of frequency adjacency. Cluster level statistics were calculated by determining the sum 
of the t-values within each cluster (cluster mass) and the maximum cluster mass 
across the data was determined. Permutations were then created by collecting the data 
from the prestimulus and task conditions across the 10 subjects in a single set. Half of 
these collected trials were placed into subset one and the remaining were placed into 
subset two. The maximum cluster mass was then determined for all possible 
permutations of the data (2 9 partitions for 10 subjects). The proportion of these 
permutations that had a smaller maximum cluster mass than the non-permuted data 
was the p-value calculated using a complete permutation test. Because of the global 
null hypothesis used for the cluster mass test, this method inherently controls for 
multiple comparisons. 

For the spontaneous phase locking analysis, two partitions were formed based on 
the ROI pair of interest and the control ROI pair comparison (i.e. for the first 
experiment, RS1-LS2 were the ROI pair of interest and RM1-LS2 was the control ROI 



pair, for the second experiment LAud-RAud were the ROI pair of interest and 
LMl-RAud were the control ROI pair). All the frequency points that were p < .05 for 
two separate criteria were determined: 1. the ROI pair of interest vs. the control 2. the 
ROI pair of interest vs. the PLV for the ROI pair of interest in noise 24 . These frequency 
points were then clustered on the basis of frequency adjacency and the rest of the 
statistics followed the procedure laid out above. 
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